%%%% mixing distributions as in Kim Shephard and Chib (1996) 
%%%% om q_j[Pr(om=j)] m_j v_j^2

ksc=[1 0.00730 -10.12999 5.79596
2 0.10556 -3.97281 2.61369
3 0.00002 -8.56686 5.17950
4 0.04395 2.77786 0.16735
5 0.34001 0.61942 0.64009
6 0.24566 1.79518 0.34023
7 0.25750 -1.08819 1.26261];